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Abstract 

Interpolating fields, used to excite hadrons from the QCD vacuum in nonper- 
turbative field-theoretic investigations of hadron properties, are explored with 
an emphasis on unconventional nucleon interpolators. The QCD continuum 
model for excited state contributions to QCD correlation functions is a central 
element in extracting the physics contained in these alternate correlation func- 
tions. The analysis confirms the independence of nucleon properties obtained 
from different interpolating fields. However, this independence comes about 
in a trivial manner. These results provide a resolution to the long standing 
debate over the optimal nucleon interpolating field to be used in QCD Sum 
Rule analyses. 

12.38.Gc, 12.38.Lg, 12.40.Yx 
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Typeset using REVT^ 



I. INTRODUCTION 



One of the instrumental operators of nonperturbative field-theoretic investigations of 
hadron structure is the hadron interpolating field. This operator is used to excite a hadron 
of specified quantum numbers from the QCD vacuum. It has long been established that there 
are two independent interpolating fields with no derivatives having the quantum numbers 
of spin 1/2 and isospin 1/2. Both are expected to excite the ground state nucleon from 
the vacuum. Various linear combinations of these interpolators are used in nonperturbative 
approaches to QCD. What distinguishes the different approaches is the manner in which the 
propagation of quarks in the QCD vacuum is determined. 

Numerical simulations of the theory via lattice regulation is the only method for probing 
deep into the nonperturbative regime of QCD. Exploitation of the Operator Product Ex- 
pansion (OPE) in QCD Sum Rules (QCD-SR) allows the near perturbative regime of QCD 
to be explored. Modeling of the QCD vacuum via instanton fluctuations in the Random 
Instanton Liquid Model (RILM) has also produced some new insights into QCD. While there 
are formal field theoretic arguments indicating nucleon properties are independent of the in- 
terpolating field, a demonstration of this in practice is an important test of these approaches 
to nonperturbative field theory. 

Some attention has been given previously to alternate nucleon interpolators on the lattice 
imj^. In these analyses conclusions were limited as the correlation functions deteriorated 
too quickly for ground state properties to be determined. In a previous paper P] it was 
established that the properties of the lowest lying state may be extracted from the first 
few points of two-point correlation functions with the use of a pole plus QCD continuum 
model inspired by QCD Sum Rule analyses. In this paper, these techniques are used to 
investigate nucleon properties obtained from correlation functions of unconventional nucleon 
interpolating fields. Some of these correlation functions suffer a loss of signal prior to a clear 
ground state domination and the QCD continuum model becomes a central element in fitting 
the correlation functions. 
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In the QCD-SR field there is a history of argument over the optimum nucleon inter- 
polating field to be used in analyses P-p]T|. This issue is recognized to be of paramount 
importance and remains unsettled ||12[. While some advocate an interpolating field for 
which the leading terms of the OPE are stationary with respect to the interpolating field 
mixing parameter [^^, others argue that a balance between OPE convergence]^ and QCD 
continuum contributions must be maintained UsHTT . 



Ideally, one would like to simply calculate with alternate interpolating fields and confirm 
that the nucleon properties remain unchanged. However, the limitations of the QCD-SR 
approach have prevented one from doing this in practice. Limitations include uncertainties in 
the values of lower dimension condensates, factorization of higher dimension operators, OPE 
truncation and convergence issues, uncertainties surrounding the role of direct instanton 
contributions to the sum rules and uncertainties in the reliability of the continuum model 
for excited states. Fortunately, the lattice approach is not plagued with the same limitations 
and the following analysis resolves this long standing debate. 

The format of this paper is as follows. Section ^ introduces the interpolating fields 
explored in this analysis, the lattice techniques used, and issues encountered in relating 
lattice and continuum (lattice spacing a — 0) formalisms. Section |ITT| highlights the QCD 
continuum model derivation. The analysis of the correlation functions is presented in Section 

for each interpolating field combination. The results are compared with other approaches 
to QCD in Section 0. Finally the conclusions regarding interpolating field invariance and 



the optimal nucleon interpolator for QCD Sum Rules are summarized in Section Vl. 



II. LATTICE CORRELATION FUNCTIONS 



^Here and in the following, "convergence" of the OPE simply means that the highest dimen- 
sion terms considered in the OPE, with their Wilson coefficients calculated to leading order in 
perturbation theory, are small relative to the leading terms of the OPE. 
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A. Interpolating Fields 



The commonly used interpolating field for the proton in lattice calculations has the form 

Xi{x) = t'^"" (u'\x)C^,d\x)) u\x) . (2.1) 



Here, we follow the notation of Sakurai C = 7472 is the charge conjugation matrix, 

a, b, c are color indices, u{x) is a u-quark field, and the superscript T denotes transpose. 
Dirac indices have been suppressed. 

In the sum rule approach, it is common to find linear combinations of this interpolating 
field and 

X2ix) = e""^ (M^«(x)Crf^(a;)) ^,u'{x) , (2.2) 

which vanishes in the nonrelativistic limit. A priori, there is no reason to exclude such an 
interpolating field [^], as the quark field operators of ( |2.2| ) annihilate the light current quarks 
of QCD. Of course, these quarks are highly relativistic when bound in the nucleon. With 
the use of the Fierz relations, the combination of the above two interpolating fields with a 
relative minus sign may be written 

Xsr(x) = 6"^"^ {u^%x)Cjy{x)) i,rd\x) 

= 2 (X2 - Xi) , (2.3) 

giving the proton interpolating field often found in sum rule calculations The alter- 

nate QCD Sum Rule interpolating field is 



Xa(x) = ^e""^ {u'\x)Ca^,u\x)) a^'^.d'ix) 
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2 (X2 + Xi) ■ (2.4) 



In this analysis we will consider both interpolating fields introduced in ( |2.1| ) and ( p.2|) 
and their interference terms such that any linear combination of these interpolating fields 
may be investigated. 



B. Correlation Functions at the Quark Level 



Hadron masses are determined through the consideration of two-point correlation func- 
tions. Here we consider the nucleon correlator 



G2{t,p) = E e-^^nr[r4 ( I T{xNix)xNm I ) 



(2.5) 



where Xn{x) may be either ( ^.1| ) or (^.2|), r4 = (1 + 74)/4 projects positive parity states for 
p = 0, and tr indicates the trace over Dirac indices. 

Correlation functions at the quark level are obtained through the standard procedure 
of contracting out time-ordered pairs of quark field operators. For the octet baryons it is 
convenient to define the correlation function 

•^('S'/i' %) = 

^afec^a'b'c'l^aa'^^^ 0) tr [Sf^ {x, 0) 5f "^(x, 0) 

^cc'T/ n\ abb' 



+S}^ (x, 0) S% ' (x, 0) S% (x, 0) 



(2.6) 



where S^'^'ix, 0) = T {g"(x), q""' (0)| and /i, /2, /s are flavor labels. For the proton interpo- 
lating field Xi of ( |2.1| ), the two-point function may be written 

G2{t,p) = e-'^-Hr [r4 r [Su, S^, CS,C-')] , (2.7) 

X 

where C = C^^. Similarly the two-point function corresponding to X2 of ( |2.2| ) may be 
written in the form 



Y-1 



(2.8) 



The interference contributions of these two interpolating fields are 



, (2.9) 
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C. Lattice Techniques 



Here we briefly summarize the lattice techniques used in the following calculations. Ad- 
ditional details may be found in Ref. 0]. Wilson's formulation is used for both the gauge 
and fermionic action (r = 1). S'f/(2)-isospin symmetry is enforced by equating the Wilson 
hopping parameters = i^d = i^- Three values of n are selected and are denoted ni = 0.152, 
K,2 = 0.154 and ^3 = 0.156. To make contact with the physical world, the mass and interpo- 
lating fleld coupling strengths calculated at the three values of k are linearly extrapolated 
to = 0.159 8(2) where an extrapolation of the squared pion mass vanishes. Differences 
between linear extrapolations to m.^ = as opposed to the physical pion mass are small and 
are neglected in the following. 

Twenty-eight quenched gauge conflgurations were generated by the Cabibbo-Marinari 
||T^ pseudo-heat-bath method on a 24x12x12x24 periodic lattice ai (3 = 5.9. Conflgurations 
were selected after 5000 thermalization sweeps from a cold start, and every 1000 sweeps 
thereafter [|T5| . 

Dirichlet boundary conditions are used for fermions in the time direction. Time slices 
are labeled from 1 to 24, with the (5-function source at t = 4. To minimize noise in the 
Green functions, the parity symmetry of the correlation functions, and the equal weighting 
of {U} and {U*} gauge conflgurations in the lattice action are exploited. The nucleon mass 
determined from xi of (|2.1| ) is used to set the lattice spacing. This estimate lies between 
other estimates based on the string tension or the p-meson mass. The lattice spacing is 
determined to be a = 0.132(4) fm and a"^ = 1.49(5) GeV. 

Statistical uncertainties in the lattice correlation functions are estimated by a single elim- 
ination jackknife |jT6|. A covariance matrix flt of the pole plus QCD continuum model over 



a range of 7 or more time slices is likely to be unreliable for 28 gauge conflgurations [17 



Instead we use standard statistical error analysis in which correlations among the flt param- 
eters are accounted for. The Gauss-Newton method is used to minimize x^- Uncertainties 
are taken from the standard error ellipse |]T8[ at = Xmin + ^■ 
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D. Operator Mixing 



The implementation of Wilson fermions on the lattice induces mixing between the com- 
posite nucleon interpolating fields of (2J.) and (|2.2| ), refiecting the breaking of chiral 
symmetry. In Ref. [1^ the mixing is argued to occur between 



0„ ^ e'^'^ (ul\x)Cd\{x)) ul{x) 
1 

(1 - 75) (Xi - X2) , 



o. 
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e'^"^ {ul\x)Cd\{x)) ulix) 
-7 (1-75) {X1 + X2) , 



(2.10a) 



(2.10b) 



and a third operator 



= e''"' (n^"(x)C7p75rf*(x)) JlYu'{x) . 



(2.10c) 



Here 



UR = ^{l+-f5)u, Ul = ^(1-75)m 



(2.11a) 



1 1 

7r = o (1 + 75) , and 7^ = - (1 - 75) . 



(2.11b) 



However, there are only two operators having isospin-1/2 and spin-1/2 and it is possible to 
demonstrate 



-20. 



(2.12) 



via Fierz transformations. For xi and X2 the expressions up to one-loop order in perturbation 
theory relating the operator matrix elements in the Pauli-Villars (PV) and Lattice (L) 



schemes are 19 



xV 



x\ 



x\ 



as 



-2 lnQV + (C{^-C2^ 
-2 lnQV + (C{^ + C2^ 



2C3L 



X2, 



(2.13) 
(2.14) 
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where 

C]^ = 37.91, = -3.21, and = -0.80, (2.15) 

for the Wilson parameter r = 1. The important point is that the interpolating field X2 does 
not mix with Xi to one-loop order. Moreover, the mixing of Xi with X2 is negligible. Hence, 
it is possible to identify the properties of these interpolating fields determined on the lattice 
with those of their continuum (a 0) counterparts to a good approximation. 

The dominant contribution to the coefficient C]" in the above expressions is from the 
self-energy corrections to the quark external lines. These corrections are accounted for in 
the mean-field improved approach, and the remaining renormalization associated with 
composite operators is relatively small. The principle renormalization constant has 
been determined in the mean-field approach and is used in the following. The nucleon 



coupling strength, Aat, is determined in absolute terms, without resorting to a ratio of the 



QCD continuum contributions as done in [pT] , p2| . In particular, the renormalization at the 
scale of 1 /a is 

Z /3k \3/2 

^Continuum Xn ^Lattice / _ " \ (2 16) 

and Z^^ = (1 — 0.73 ay) ~ 0.80 at (3 = 5.9. The k dependence of this wave function 
renormalization is very different from the naive normalization 

^Continuum _ ^ ^Lattice ^2^)^/^ , (2.17) 

a^/2 

and is crucial to recovering the correct mass independence of the Wilson coefficient of the 
identity operator. 

III. THE QCD CONTINUUM MODEL 

Here we briefiy review the QCD continuum model implementation in Euclidean space as 
examined in detail in Ref. We start with the two-point correlation function of ( |2.5| ). At 
the phenomenological level, one inserts a complete set of states A^* and defines 
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{0\xN{0)\N\p,s)=y^u{p,s), (3.1) 

where the couphng strength, A^, measures the abihty of the interpolating field xn to an- 
nihilate the i'th excitation with nucleon quantum numbers. For p = and Euclidean time 
t ^ oo, the ground state dominates and G2{t) A^e^^^^*. The spectral representation is 
defined by 

POO 

G2{t) = / p{s)e-''ds, (3.2) 

JO 

and the spectral density is, 

p{s) = Xl6{s-Mr,) + as) (3.3) 

where ({s) provides the excited state contributions. 

The form of the spectral density used in the QCD continuum model is determined by 
the leading terms of the OPE surviving in the limit t — > 0. Here, the combination of 
interpolators XiXi is considered. The derivation of the QCD continuum model contributions 
to other correlation functions proceeds in an analogous fashion. In Euclidean space, G2{t) 
has the following OPE 

^ , , 3 ■ 5 V 1 28 m„a 14 m^a^ 

-^i^ + --). (3.4) 

The spectral density used in the QCD continuum model is defined by equating ( |3.2| ) and 
( p.4|) . The QCD continuum model is defined through the introduction of a threshold which 
marks the effective onset of excited states in the spectral density. Keeping the first two 
terms of ( p.4|) , the phenomenology of 6*2 (^) is 



roo 

G2it) = Xje'^^^' + ^ p{s)e-''ds, (3.5a) 

J so 

= A?e-*^-* + e|^e--*x (3.5b) 

(f 1 , ^0 , 1 ^0 , 1 si 1 1 sl\ 28m,a f 1 sp 1 si 1 si 1 s^]] 

[\t^ 2 Q 2A 120 t \ 25 \ 2 Q 2A t \ ' 



The parameter ^ governs the strength of the QCD continuum modeL In the continuum hmit 
{a —>■ 0) ^ = 1 but here is optimized with X^, Mjv, and Sq to account for enhancement of 
the correlator in the short time regime reflecting lattice anisotropy. C, is an overall QCD 
continuum model strength and is expected to be independent of the quark mass. With 
this approach, the effects of lattice anisotropy may be absorbed through a combination of a 
larger QCD continuum model strength > 1) and marginally larger threshold (sq). 

Infrared lattice artifacts are not a significant problem for this approach as the Fourier 
transform weight exp{—ip-x) is correct for all propagator paths including those which wrap 
around the lattice spatial dimensions. The ultraviolet lattice cutoff may be modeled in a 
manner similar to that for the QCD continuum model. However the modeling becomes 
insignificant by the second time slice following the source. Instead we simply discard the 
source and first time slice when fitting the correlation functions. 

IV. LATTICE CORRELATOR FITS 
A. ( XlXi ) Correlation Function 

Xi is the standard nucleon interpolating field used in lattice analyses. Its overlap with 
the nucleon ground state is excellent. This lattice correlation function is fit with ( |3.5| ) in 
a four parameter search of X^, M^, sq and ^ in analysis intervals from t = 6 — tj where 
tf ranges from 11 through 23. Figure || illustrates the lattice data and these 13 fits at our 
intermediate value of quark mass. Similar results are seen for k = 0.152 and 0.156. 

The similarity of the 13 pole plus QCD continuum fits establishes that the QCD contin- 
uum model effectively accounts for excited state contaminations in the correlation functions 
and allows the extraction of the ground state properties from a regime as small as t = 6 ^ 11. 
For an in-depth examination of this correlator see Ref . . The quark mass dependence of A^v 
is illustrated in Figure 0. Table | summarizes the fit parameters for the regime t = 6 — > 20. 
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B. ( X2X2 ) Correlation Function 



Unlike QCD Sum Rule analyses, correlators may be determined on the lattice for any 
interpolating field without regard to the restrictions of OPE convergence issues or operator 
factorization assumptions. In constructing the QCD continuum model, only the first few 
terms of the OPE are required. The OPE for the interpolating fields X2X2 

3-5^/1 4 m„a 2 mW 



G2{t) 



2%4 5 t5 5 ^4 

87r^ ( qq ) 
15 t3 

and the QCD continuum model is derived in an analogous manner to that outlined in Section 



+ ^<M^ + ..A. (4.1) 



PTT] . The phenomeno logical side of G2{t) is 



' 1 1 1 1 I 1 gp] 4mga fl so ls| ls| ls|' 

t"^ 2 t'' Q 2A 120 t \ 5 2 Q 24: t 



Figure |^ illustrates the lattice correlation function and the final fit. The choice of r4 in 
( [^.5D projects out positive parity nucleon states when p = and therefore the correlation 
function must remain positive. At t = 13 the lattice correlation function data changes sign, 
and indicates a loss of signal. 

The fit from t = 6 ^ 12 using a pole plus QCD continuum model leads to fit parameters 
where the pole lies above the continuum threshold. The position of the pole is insignificant 
as its removal has little effect on the x^/dof. Fixing the pole at the previously determined 
nucleon masses returns an optimum value for \\ which is negative, and once again unphysical. 
The fit illustrated in Figure |^ employs the QCD continuum model alone. Hence there is 
no evidence of any overlap of X2 with the ground state nucleon in this correlation function. 
While the results illustrated here are for our intermediate value of quark mass considered 
on the lattice, similar results are seen for the lighter and heavier quark masses. The fit 
parameters are summarized in Table y. The QCD continuum threshold is not too different 
from that for XiXi- 
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Figure ^ illustrates the quark mass dependence of ^. As anticipated, ^ is independent of 
niq. This quark mass independence confirms the negative sign of the rriq correction appearing 
in the OPE for X2X2 ^"^^ ^se of mean-field improved operators. It also confirms the 
perturbative role of the quark mass operator. While the heaviest current quark mass used 
in this investigation is similar to that of the strange quark, it is still light on the scale set 
by the nucleon mass. 

It is also interesting to note that ^ is much closer to 1 than for the correlators of XiXi- 
That this might be the case is eluded to by the opposite signs of the leading terms of the 
OPE in (|4.1|) , providing the possibility of cancelations in the short-time perturbative regime 
of the correlation function. 



C. ( X1X2 + X2X1 ) Correlation Function 

Since the square of A2 is small, one might be able to recover a signal for the overlap of 
the nucleon ground state and X2 by considering the correlation function for XiX2- Figure |^ 
illustrates a three parameter fit from t = 6 ^ 12 using the pole plus QCD continuum model 
derived from the OPE for ^{xiX2 + X2X1) 



G2{t) 



3-5/1 4 rriga 2 m^a^ 



Stt^ (qq) 



In this fit, the nucleon ground state pole position has been fixed at the previously deter- 
mined nucleon masses. The fit parameters are summarized in Table While there is 
sufficient information in the correlation function to determine a value for the nucleon mass, 
the corresponding uncertainties are large. 

The leading terms of the OPEs for X2X2 |(xiX2 + X2X1) given in (|4l| ) and (|4.3| ) are 



equivalent up to a normalization factor of five. Since the continuum model is constructed 
to accommodate these leading terms, one expects the QCD continuum model parameters. 
So and ^, for these two correlators to be similar. A comparison of Tables |I| and |II| indicates 
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that this is indeed the case. 

Figure ^ illustrates the expected quark mass independence of ^, again confirming the 
negative sign of the rriq correction appearing in the OPE for |(xiX2 + X2X1) the use of 
mean-field improved operators. Similarly ^ ~ 1 as anticipated by the opposite signs of the 
leading terms of the OPE of p^). 

The linear extrapolation of (AiA2)^^^ to is illustrated in Figure ^ The combination 

1 /2 

(A1A2) shows little sensitivity to the quark mass. This contrasts the dependence of Ai il- 
lustrated in Figure H, where Ai decreases as the quarks become lighter. Thus, A2 increases for 
decreasing quark mass. This reflects the fact that X2 vanishes in a nonrelativistic reduction. 
At the chiral limit (AiA2)^^^ = 0.0014(10) GeV^. Systematic uncertainty in the extrapo- 

1 /2 

lated value of (A1A2) may be estimated using the quark mass dependence suggested by 
Chiral Perturbation Theory p3[ as in 0. Here the systematic uncertainty in extrapolating 



is negligible relative to the statistical uncertainties. 

With the previous result Ai = 0.013(2) GeV^, the nucleon coupling strength for X2 is 
found to be A2 = 0.00016(22) GeV^, approximately 100 times smaller than Ai. In short, 
there is only one nucleon interpolating field that has significant overlap with the nucleon 
ground state, namely 

Xi(x) = e'^''^ (u'^''{x)C^,d\x)) u%x) , (4.4) 

and 

AsR^ Aa ^ 2Ai = 0.026(4) GeV^ , (4.5) 
A2^0. (4.6) 



The results for all the considered interpolating fields are summarized in Table IV |24 . 



V. COMPARISON WITH OTHER CALCULATIONS 



It is worth noting that the nonperturbative QCD Sum Rule predictions for Asr have 
remained quite stable over the years, despite the fact that the early calculations have a 

13 



number of numerical errors in Wilson coefficients and anomalous dimensions |^5|. This 
cannot be said for the more model dependent predictions. It is ironic that, in some cases, 
the model calculations were pursued due to reservations about the reliability or validity of 
the QCD Sum Rule approach. Table |V| summarizes a collection of predictions taken and 
updated from Ref. |]2B| and 



QCD Sum Rule predictions of A2 or Aa are more uncertain. This is largely due to a lack 
of rigor in the analysis of the sum rules. Often, the region of validity in Borel space is simply 
postulated with little regard to OPE convergence or the size of continuum model contribu- 
tions. Many authors have fixed the continuum threshold to a preferred value or excitation 
energy rather than leaving it as a search parameter to be optimized. The upper limit of the 
Borel region must be monitored as it is a function of the three required fit parameters, M^v, 
Atv and Sq and varies for different interpolating fields. The failure to monitor these issues in 
existing analyses, is largely responsible for the apparent inconsistencies between sum rules 
derived from different interpolating fields. 

Ref. [^] is one of the few sum rule analyses where these issues are rigorously implemented. 
However, the interpolating field X2 was not considered there. The sign of the quark con- 
densate term in ( |4.1| ) indicates the two sum rules will be saturated by both positive and 
negative parity states. A careful analysis of these sum rules has not yet been attempted. 

In the QCD-SR discussion of Ref. it was concluded that the overlap of xa with the 
nucleon, Aa, must be negligible, due to the vanishing of most of the Wilson coefficients to 
dimension 8. However, this conclusion need not be the case. Higher order terms of the OPE 
starting at dimension 9 are not zero and could easily give rise to large overlap with the 
nucleon as discovered here. The only conclusion that may be drawn from these sum rules is 
that the pole contribution is small relative to the QCD continuum model contribution. The 
pole contribution is not necessarily small in absolute terms. 

The first five entries of Table [V] summarize results for Asr obtained from the consideration 
of two-point correlation functions, and these compare favorably. The same cannot be said 
for Aa- The RILM prediction |^ is Aa = 0.040(2) GeV^ and is large compared to the lattice 
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prediction of Aa = 0.027(5) GeV^. Figure 3 of Ref. displays significant discrepancies 
between a global fit of the six nucleon correlators considered and two of the correlators. 
These two correlators are both dependent on Aa and these discrepancies are not reflected in 
their quoted uncertainty of ±0.002 GeV^. Their conclusion "What is even more important, 
the simple 'nucleon pole plus continuum' model gives a very good simultaneous description 
for the complete set of correlation functions" is difficult to justify in the RILM, particularly 
in light of these new results. 



VI. CONCLUSIONS 

A. Interpolating Field Invar iance 

Ground state nucleon properties are independent of the interpolating field used to ex- 
cite the baryon from the vacuum. This invariance is satisfied in a trivial manner. The 
interpolating field which vanishes in the nonrelativistic limit, 

X2{x) = e""^ (u'\x)Cd\x)) 75M^(x) , (6.1) 

has negligible overlap with the nucleon ground state. Inclusion of X2 components in interpo- 
lating fields only increases the statistical uncertainties of lattice QCD correlation functions. 



B. Optimal Interpolator for QCD Sum Rules 

This analysis indicates that, to a good approximation, X2 excites pure QCD continuum. 
Since X2 has negligible overlap with the ground state nucleon, it is tempting to simply 
conclude that the optimum interpolating field is xi- While this is certainly the case for 
lattice QCD investigations, it is not obviously the case for QCD Sum Rule analyses. 

The optimal nucleon interpolator must involve Xi ^is this interpolator is required to 
maintain overlap with the ground state. The task is to determine the optimal mixing of X2- 
The Borel improved QCD Sum Rules for the generalized interpolator 
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Xo = Xi + (3X2 , 



(6.2) 



are 



where 



5 + + ^6 ^-4/9 

64 

5 + 2^3 + 5^ 

256 



7 - 2/3 - 5p'' 2.4/9 13-2/3- 11/32 mgg^ 
^ 24 96 



(6.3a) 



and 



7 - 2/3 - 5/32 
16 

3(1 -/32) 



16 

3 + 2/3 - 5/32 



1 _ e-^oA^' 



a b 



Xl Mm e-*^^/^-^ - Xl, Mm* e"^^*/^^ 



0, IVlN* 



(6.3b) 



-(27r)2 ( gg )= 0.450 GeV^ 



a 



b = (27r)2 ( Gl^C^"" )= 0.474 GeV 



TT 

mg = -114^:^ = 0.65, 
^ _ log(M/AQCD) 

log(/i/AQCD) 

Ao = (27r)^Ac) . 



(6.4a) 
(6.4b) 
(6.4c) 

(6.4d) 
(6.4e) 



Here M is the Borel parameter and plays a role similar to the inverse Euclidean time of 
the Lattice approach. The condensate values are taken from Ref. where = 0.5 GeV 
and Aqcd = 0.1 GeV. The continuum model contributions are indicated on the left hand 
side of the sum rules where they appear in brackets as subtractions from the terms of the 
OPE surviving in the limit M 00. To aid the following discussion, both a positive parity 
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ground state and a negative parity excited state are included on the right-hand-side of the 
sum rules. 

The first QCD-SR of ( |6.3a| ) is known to have uncontrollably large perturbative corrections 



to the Wilson coefficient of the identity operator [^. In leading order, these corrections are 
independent of (3 and are approximately 50%. As a result this sum rule must be discarded 



In the QCD-SR approach, approximations are made at both the quark level and the phe- 
nomenological level. At the quark level, the OPE is truncated and the Wilson coefficients 
are calculated perturbatively. This sets a lower limit for the Borel mass. At the phenomeno- 
logical level, the spectral density is approximated by a pole plus the QCD continuum model. 
Maintaining ground state dominance on the phenomenological side of the sum rule sets an 
upper limit on the Borel mass. By including X2 components in an interpolating field, one 
can reduce the continuum contributions excited by xi allow a broader Borel analysis 
window. 

One of the most difficult things to monitor in the QCD-SR approach is whether the 
OPE is sufficiently convergent for a particular value of Borel mass. The lattice results 
presented here indicate the X2X2 correlator has the fastest converging OPE, as its overlap 
with the nucleon ground state is negligible. Similarly, the combination XiXi produces an 
OPE with the slowest convergence, as this correlation function is dominated by the ground 
state nucleon for small Borel masses. 

Hence, errors made in truncating the OPE are dominated by errors in the XiXi com- 
ponent of the general correlator. The relative error in the OPE truncation can be reduced 
by adding X2 components to the correlator. However, the X2 components in the OPE are 
simply subtracted off again by the continuum model terms. Hence the relevant error is the 
absolute error. For \j3\ < 1, this error is dominated by XiXi components of the correlator. 
As a result, OPE truncation errors are approximately independent of (3. This crucial point 
has been neglected in previous arguments regarding the optimal nucleon interpolating field. 

Since X2 has negligible overlap with the nucleon, the ground state contribution is also 
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independent of (5. Hence, the size of the continuum model contributions is the predominant 
factor in determining the optimal interpolator. Figure |^ illustrates the contributions of the 
continuum model terms in ( |6.3b| ) for M = 0.938 GeV and sq = 1.4 GeV. The following 
discussion is not dependent on the precise values of these parameters. The first point to 
be made is that contributions from the continuum model are largest for f3 ~ —0.2. This 
selection of mixing is the worst possible choice for extracting information on the ground 
state nucleon. 

Figure || also indicates it is possible to have vanishing continuum model contributions at 
f3 ~ —1.5 or (3 = 1. However, we are relying on the continuum model to account for strength 
in the correlator that does not have its origin in the ground state. Without a continuum 
model, one would need to include additional poles on the right-hand-side of ( |6.3b|) to account 
for positive and negative parity excitation strength. For /3 < — 1.5 or /3 > 1.0 the correlator 
is negative indicating the sum rule is saturated by a negative parity state. 

Thus the optimal interpolator is /3 ~ —1.2 or /9 ~ 0.8. To discriminate between these 
two regimes, we turn to the higher dimension operators (HDO) which do not contribute to 
the continuum model. It is these terms that provide crucial information on whether the 
strength in the correlator lies in the ground state or the excited states. If these terms are 
absent, the optimal fit of the correlator is obtained when Ao ^ and sq — > 0. In this case 
the continuum model becomes the Laplace transform and the fit is perfect. Hence the HDO 
terms should be large in magnitude. A change in sign from the leading terms of the OPE 
will also assist in distinguishing ground state strength from excited state strength as the 
change in the curvature of the correlator will be more prominent. 



The last term of ( |6.3b|) is a HDO term, and its value is plotted as a function of beta in 
Figure The HDO contributions are larger for j3 ~ —1.2 than for j3 ~ 0.8. In addition, 
the sign of the HDO contribution is opposite that of the continuum model contributions. 
Hence the preferred regime is /3 ~ —1.2. In fact, optimization of the three fit parameters, 
Mtv, Xo and Sq {\o* = 0) of ( |6.3b| ) for /3 = 0.6 ^ 0.8 results in fit parameters describing 
pure continuum with sq ~ M^- Information to separate the ground state pole from the 



continuum is insufficient for (3 = 0.6 — * 0.8. 

In summary, the lattice results indicate OPE convergence, and ground state pole contri- 
butions are approximately independent of (3. Considerations of the size of continuum model 
contributions and the sign and magnitude of HDO operators leads to the preferred value of 

/3 = -1.2 ±0.1. (6.5) 

A more precise determination of (3 will depend on the details of limits for continuum model 
contributions, HDO values, condensate values, and other parameters of the sum rules. 
Hence, this analysis supports the selection of /3 = — 1 |l7HTI|, over (3 = —0.2 Q-Q. At 



P = —0.2, where the leading terms of the OPE are stationary with respect to P the 
continuum contributions are maximal. The positive value and small magnitude of the HDO 
indicates the stability of the leading terms of the OPE will not be realized as stability in 
the ground state mass, coupling, nor in the continuum threshold. 



C. Future Investigations 

These techniques may be used to determine the optimal interpolating field for any sum 
rule involving Xi and X2 components. Each sum rule will have an optimal selection for f3. 
The overlap of spin-1/2 and spin-3/2 interpolating fields are known to yield nucleon sum 
rules which offer stability in the fit parameters that cannot be obtained from the more 
common sum rules considered here. It will be interesting to discover if the historical selection 
of /3 = — 1 is indeed optimal. 

While it is important to establish the optimal mixing of interpolating fields for QCD-SR 
analyses, one should not overlook the fact that there is a range of values for (3 where the 
sum rules are expected to work. Moreover, the ground state contribution to all these sum 
rules is equivalent to the 1% level. In other words, the right-hand-side of ( |6.3b| ) for a single 
pole plus continuum model is independent of j3. After the first sum rule is written down, 
additional sum rules may be introduced with merely one new fit parameter (sq) per sum 
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rule. Since direct instanton contributions to the sum rules are not independent of (3 []3T| , 



one has an excellent opportunity to see if direct instanton contributions really are necessary 



to maintain sum-rule consistency p^] 



Future lattice QCD investigations should aim to make a direct comparison of the OPE 
and lattice correlation functions. The Wilson coefficients and vacuum expectation values of 
normal ordered operators could be determined directly from OPE fits to the lattice data. 
Such a comparison would test the validity of the OPE in the nonperturbative sector and our 
understanding of quantum field theory |l33| . 

A direct comparison of lattice and continuum formalisms requires the use of an im- 
proved or perfect |^ lattice action to reduce or eliminate lattice anisotropy in the short 



time regime of lattice correlation functions. Alternatively, the Wilson coefficients of the 
Euclidean-space correlation function may be derived via lattice perturbation theory. 

An extremely fine lattice spacing is required to provide a sufficient number of lattice sites 
within the radius of convergence of the OPE. In the most optimistic case, the two invariant 
nucleon sum rules of a given interpolator could be isolated such that, to dimension eight, each 
correlation function would have up to four parameters to be determined when extracting 
OPE coefficients. Ultraviolet cutoff considerations, OPE convergence issues and the need 
for error estimates in the fit parameters place the lattice spacing at less than 0.05 fm. 

Essential information on the importance of direct instanton contributions to the OPE 
can be obtained from such an investigation. The OPE coefficients and vacuum expectation 
values are determined first by matching OPE and lattice correlation functions in which 
there are no direct instanton contributions. This approach determines the OPE in a self 
consistent manner. Then, other correlators in which direct instanton contributions are 



argued to be important ||3T| can be examined. Discrepancies between the OPE and lattice 
correlation functions would signal the possible importance of direct instanton contributions. 
If existing predictions for direct instanton contributions to the OPE resolve the differences 
in the correlators, then one has compelling evidence of a non-trivial role for direct instanton 
contributions to the OPE. The importance of such investigations warrants further effort in 
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this direction. 
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TABLES 



TABLE I. ( XiXi ): Four parameter search for the pole plus QCD continuum model. 



Parameter 


Kl — U.iOZ 


K,2 — U.i04 


Av3 — U.iOO 


^cr — U.loy o[Z ) 


Mjva 


1.109(8) 


0.983(8) 


0.858(8) 


0.628(17)^ 


Aia^ (xlO-2) 


1.17(5) 


0.94(4) 


0.75(3) 


0.38(7) 


Sod 


1.68(3) 


1.58(3) 


1.49(4) 


1.32(7) 




6.83(10) 


6.74(9) 


6.62(9) 


6.42(19) 


^ from OPE fit 


5.3(1) 


5.6(1) 


5.8(1) 




^The physical prot 


on mass sets the lattice 


spacing a = 0.132(4) fm. 




TABLE IL ( X2X2 )• Two parameter search for the pure QCD continuum model. 


Parameter 


Ki = 0.152 


K2 = 0.154 


K3 = 0.156 


Kcr = 0.1598(2) 


soa 


1.58(2) 


1.48(2) 


1.40(2) 


1.23(5) 




1.55(4) 


1.55(4) 


1.56(4) 


1.56(8) 


^ from OPE fit 


1.36(6) 


1.34(6) 


1.35(6) 




TABLE III. i 


( X1X2 + X2X1 )'■ Three parameter search for the pole 


plus QCD continuum 


model. Mjv has been fixed to the previously determined lattice values. 




Parameter 


Ki = 0.152 


K2 = 0.154 


K3 = 0.156 


Kcr = 0.1598(2) 


(AiA2)'/'a3 (xlO- 


3) 0.50(14) 


0.42(15) 


0.47(14) 


0.43(30) 


soa 


1.55(5) 


1.45(6) 


1.36(7) 


1.18(14) 




1.57(4) 


1.57(4) 


1.57(4) 


1.57(10) 


$, from OPE fit 


1.42(4) 


1.41(5) 


1.42(4) 
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TABLE IV. Summary of Lattice Results for the pole plus QCD continuum model. 



Interpolating Fields 






K 


So 










(GeV) 


= e™<^ \u"' C 75 a 1 u 


0.938'^ 


r\ r\-\ /f^\ 

0.013(2) 


6.42(19) 


1.98(11) 


_ abc ( a fj jb\ c 




nom 


1 5fil'7'> 




\ {XiX2 + X2X1) 


Fixed 


0.0014(10) 


1.58(9) 


1.76(20) 


XSR = (n'^C7^n'') 757'^fi'^ 


0.96(3) 


0.027(5) 


4.61(14) 


1.92(11) 




0.91(3) 


0.022(5) 


3.54(10) 


1.82(10) 



^Defines the lattice spacing a. 

'^Inferred from XiXi ^^id ^ {X1X2 + X2X1) results. 



TABLE V. Comparison of predictions for Asr for various approaches to QCD. 



Approach 


Reference 








(xlO-2 GeV3) 



Lattice (mean-field improved) 
Lattice (conventional renormalization) 
Lattice (coordinate space) 
QCD Sum Rule 
Instanton Liquid 
Baryon wave functions {pc^ 
Quark Model 
Bethe-Salpeter amplitude 
Quark Model 
MIT Bag Model 



0) 



This work 



Gavela et al. |27| 



Chu et al. |22| 
Leinweber |p 
Schafer et al. 
Brodsky et al. ]2( 
Thomas and McKellar [| 
Tomozawa [ p^ , |39| 
Milosevic et al. [|0) 
Donoghue and Golowich 



2.7(5) 
2.4 
2.2(4) 
3.1(6) 
3.2(1) 
12 
8 

2.5 
2 

1.27 
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FIGURES 

FIG. 1. The two-point correlator at k = 0.154 for the nucleon interpolating fields XiXi of 



(2.1). The fits for the 13 analysis intervals are illustrated. The source position is at to = 4. Neither 



the source nor t = 5 are included in the fit. 

FIG. 2. The quark mass dependence of the nucleon coupling strength Ai corresponding to the 
interpolating field of (p.ll). 



FIG. 3. Lattice correlation function for interpolating fields X2X2- '^^^ illustrated fit employs 
the QCD continuum model alone. 

FIG. 4. The dependence of ^ on the quark mass. The displayed independence confirms the 



sign and magnitude for the Wilson Coefficient of the mq term in the OPE of (4J.) for X2X2- 



FIG. 5. Lattice correlation function for interpolating fields ^(xiX2 + X2Xi)- The nucleon mass 
has been fixed at previously determined values for this fit. 

FIG. 6. The dependence of ^ on the quark mass. The displayed independence confirms the sign 
and magnitude for the Wilson Coefficient of the ruq term in the OPE of (|4.3| ) for ^{xiX2 + X2Xi)- 



FIG. 7. Linear extrapolation of the coupling strength (AiA2)^^^ to Kcr- The y-axis scale is one 
tenth of that in Figure ^. 

FIG. 8. Continuum model (solid curve) and higher dimension operator (HDO) (dashed curve) 
contributions to the Borel improved QCD-SR of ( 6.3b| ) plotted as a function of the interpolating 
field mixing parameter f3. 
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